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jJIj We report measurements of mesoscopic fluctuations of elastic cotunneling in Coulomb 

\Q blockaded quantum dots. Unlike resonant tunneling on Coulomb peaks, cotunneling in the 

^ valleys is sensitive to charging effects. We observe a larger magnetic field scale for the 

^ cotunneling (valley) fluctuations compared to the peaks, as well as an absence of "weak 

localization" (reduced conductance at B=0) in valleys. Cotunneling fluctuations remain 
correlated over several valleys while peak conductance correlation decreases quickly. 

o 

g The Coulomb blockade (CB) in small metalUc and semiconductor structures provides a system in 

I which transport is dominated by the effects of electron-electron interactions [1, 2]. At low temperature, 

^ transport in small systems is also influenced by quantum interference, leading to mesoscopic fluctuations 



with universal statistical features [3, 4]. Recent theoretical work has provided a detailed understanding of 
mesoscopic fluctuations for the case of non-interacting electrons and has called attention to the connection 
between the universal statistics and quantum chaos [5-8]. Coulomb blockaded quantum dots provide a 
useful system to move beyond the non-interacting picture, yet remaining relatively simple: interactions can 
often be treated in terms of a single charging energy. 

The dominant experimental signature of the CB is the appearance of narrow peaks in conductance 
as an external gate voltage is swept, changing the potential of the dot. Conduction is suppressed between 
peaks when the temperature, T, and voltage bias, V^^, are less than the energy = e^/c required to add 
one electron to the dot. In the quantum regime, (kT, Vgjj) < A « , (where A = iTth^ / m*A is the mean 
level spacing and A is the dot area) transport on CB peaks is mediated by resonant tuimehng and is insen- 
sitive to charging effects. This has been demonstrated, for instance, in recent measurements of peak 
height statistics [9, 10] which gave excellent agreement with single-particle random matrix theory (RMT) 
[11-13]. In contrast, conduction between CB peaks is mediated by so-called cotunneling, which is sen- 
sitive to charging [14-17]. At low temperatures, the dominant cotunneling process is elastic and, like the 
peaks, exliibits interference effects in the form of random but repeatable conductance fluctuations [18]. 
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In this Letter, we report the first measurements of mesoscopic fluctuations of elastic cotunneling 
in the valleys between CB peaks. The quantum dots used are ballistic, with a "chaotic" shape allowing 
comparison to a universal theory [18] and make use of shape-distorting electrostatic gates (Fig. 2(b), in- 
set) to allow ensemble statistics to be gathered on a single device. We observe an increase in the 
characteristic magnetic field scale of conductance fluctuations in the valleys compared to peaks in the CB 
regime and find that valleys do not exhibit the decrease in average conductance at B=0 ("weak localiza- 
tion") that is seen in peaks. We also compare cross-correlations of neighboring peaks and valleys, finding 
longer correlation among valleys. This effect is explained in terms of the number of quantum levels par- 
ticipating in transport. 

In the valleys between CB peaks, conductance is mediated by two mechanisms: elastic and inelas- 
tic cotunnehng [16] . Each contributes to the average conductance, 

(G.,) = ^^f- + -1, (G,,) = ^^(ytr)<^+^T, (1) 

where £'f,(/,) is the difference between the Fermi energy in the leads and the next available state for elec- 
tron (hole) tunneling (i.e. Eg + E^ = E^, and E^ = Ef^ = E^/2 in the middle of the valley, see Fig. 2(c), 
inset). At low temperatures, kT < -^JE^A , conductance is dominated by the elastic mechanism, consisting 
of virtual tunneling of charge over an energy barrier of height ~E, the smaller of {E^,E^). Near the val- 
ley center, where E ~ E^I2 » A, virtual tunneling takes place through a large number -E/A of levels. 
Conduction through each level fluctuates randomly with external parameters, resulting in mesoscopic co- 
tunnehng fluctuations that do not average away even for large A [18]. 

Experimentally, one can distinguish resonant tunneling fluctuations (on peaks) from cotunneling 
fluctuations (in valleys) by the characteristic magnetic field scale, B^, defined as the width of the autocor- 
relation Q^iiAB), where Qj(AB) = (^,(5)^y(fi + A5))^/(Vvar|,-^var|^], and ^ = G-(G)g. In 
particular, 5™"^^ is the field required to pass roughly one flux quantum, (^^ = hie, through a typical area 
difference accumulated by chaotic trajectories in the time T ~ hjE limited by virtual tunnehng at energy E, 
giving 

Bj-'''^y = K{^JA)4EIE,, (2) 

where Ex = hvp/A^'^ is the ballistic Thouless energy and ?c is a device-dependent geometrical factor 
[18]. Equation (2) applies when most virtual ttajectories are fully chaotic, E < E^. In the opposite case, 
E > Et, Ej replaces E, giving Bj"^^^ ~K(I>^/ A. For comparison, the characteristic field scales of con- 
ductance fluctuations on CB peaks and in open quantum dots are 

S/'^* = K{^, IA)4KlE, , B;"^" = K{^JA)^Y,^,IE,, (3) 

where Ftot = Fi + Fr + is the total broadening due both to escape {Fi,. = {hAl7te^)Gir) and 
dephasing [7, 19, 20]. One expects (B/'"'' , B^"P''") < B/""^ from Eqs. (2) and (3) since 
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{A,r,gi) < E ~ E^ll. Equation (2) further implies that 5^™"^>' will have a maximum at mid-valley, where 
E = EJ2, and decrease to match B/'^"'^ as the peak is approached and £ — > . Making the approxima- 
tion (G^,)oc(£^-i + £;j-i)oc~ £-1 and noting that B/""'^ oc ^ further unpUes that B^^{G^i) ~ 
constant. 

We report measurements for three quantum dots defined using Cr/Au electrostatic gates on 
GaAs/AlGaAs heterostructures. Important parameters for the dots are given in Fig. 4(d). The dots use 
similar designs with adjustable point contacts and two or three shape-distorting gates as shown in Fig. 
2(c). Measurements were made in a dilution refrigerator using an ac voltage bias with V^^ ~ 5 |J.V (< A, 
kT) at 13Hz. The electron temperature, T~ 100 mK, measured from CB peak widths, easily satisfies the 
requirement kT < -JE^A ~ 600 mK. Conductance fluctuation statistics were measured by rastering over 
gate voltage, V„, and field, B, to yield a 2D scan of conductance, as seen in Fig. 1(a), then finding the 
maxima of peaks from cosh'^ fits [21] and the minima of valleys from parabolic fits. Random conduc- 
tance fluctuations are seen in the peaks, valleys, and all interpolated paths (see Fig. 1(b)). The field scale 
of fluctuations is defined by the condition C,- ,(0.325 B^.) = 0.82. This somewhat awkward definition was 
chosen to coincide with Ref. [18] while allowing B^ to be measured in the universal region, 
Ci i(AB ~ 0) ~ 1, where presumably nonuniversal features of the dot geometry such as short trajectories 
are not important. The periodic change in B^ from a maximum in the valley to a minimum on peak can be 
seen in the width of the Q iiAB) color scale. Fig. 1(c), or directly from the autocorrelations shown in 
Fig. 1(d). We note that the expected difference in ihs functional forms of C, ,(iiB) for peaks versus val- 
leys [12, 13, 18] is not resolvable in the data. 

The ensemble-averaged field scale, B^ , and conductance, {G)g , were obtained from statistically 
independent peak-valley-peak data sets, sampled by changing the voltage applied to one of the shape- 
distorting gates. Figure 2(a) shows B^ and (G)^ for 31 peaks and 22 valleys from dot 1, estimated to 
contain -14 statistically independent data sets. The measured values, B^ ^"'^^'^y_ ^ ^ + 0.3 mT and B^ 
= 4.0 ± 0.2 mT, give a ratio B^ P'"^'' _ j consistent with ratios of 1.4 - 1.8 for the other dots. 

While the general behavior of B^ in valleys compared to peaks agrees with theoretical expectations, the 
ratio is smaller than expected from Eqs. (2) and (3). In particular, one would expect ^J'^^^^yjg^P^"^ ^ 
(18mT /4.6 mT) ~ 4 for dot 1. This discrepancy is found in all of the dots and is not understood at pres- 
ent. Figure 2(b) shows the approximately linear dependence of B^. on {G)g for the data in Fig. 2(a). 
We note, however, that this relation assumes only a single carrier — either electrons or holes — and so can 
only be applied near the valley bottom or peak top. 

To show that 5™"^^ is enhanced due to charging effects, it is useful to compare B^. for both 
peaks and valleys to the characteristic field B^"''^" for open dots (G^^f > e^/h). Ensemble-averaged field 
scales, B^"'^™ (using the definitions of CijiAB) and 5^ above), obtained from -30 statistically inde- 
pendent conductance traces at three different lead conductances in dot 1 are shown in Fig. 2(c) along with 
^^peak ^ ^2^^^ ^jaiiey ^ ^2/^^ tunneling regime. These data show 

that Bp "^^^ converges to B^.^^"^ as the dot becomes isolated, while g^^"''^^ is considerably larger than 
either. This supports the observation that the characteristic (single-particle) energy scales for transport 
though open dots and for resonant tuimeling, T and A, coincide at the onset of blockade, and that both are 
smaller than E (set by classical charging) which determines B^. ^"^^^^ . 
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It is known theoretically [11, 22] and from recent experiments [9, 10] that the normalized on-peak 
conductance, g{B) = G{B) j {G) , is lower when time-reversal symmetry is obeyed (i.e. at B=0), in 
analogy to weak localization in open dots [23, 24] and ID and 2D disordered conductors [25]. In con- 
trast, elastic cotunneling in the valleys is not expected to show weak localization — that is, g^^n^y should 
not be suppressed at B=0 [18]. We have investigated the change in average conduc- 
tance, 8g = {g{B > Bf.) - g(B = 0)) , for 8 1 independent pairs of peaks and valleys measured in dot 3 . As 
shown in Fig. 3, we find Sgp^ak = ^-^'^ for the peaks, somewhat smaller than the RMT re- 
sult, Sgp^ak = 1/4 [11, 22]. Since the RMT calculation assumes F « kT « A while the measurement 
has F ~ 0.7A and kT ~ A, it is reasonable that theory should overestimate the measured value. For the 
valleys, we find that gy^Uey lacks any significant dip on a field scale of B^ (~ 8 mT for dot 3) around B=0 
(Fig. 3(b)), in agreement with theory [18]. We note that averages such as g{B) are difficult to measure in 
the CB regime because, unhke in open dots, fluctuations in g are on the order of g itself [11, 18, 22], as 
illustrated in Fig. 3. 

Finally, we investigate correlations between neighboring peaks and valleys as a function of sepa- 
ration (in units of peak spacing), AA^. Ensemble-averaged cross-correlations, C(AB, AA') = Qj+^iAB), 
for 5 peaks and valleys are shown in Fig. 4. The maximum of C(AB, AN) at AB = is seen to decrease 
to ~ for AA'^ > 2 for the peaks, whereas for valleys the correlation remains high, C(0,3)~ 0.5. This is 
also seen in Fig. 4(c) for C(0, AA'^) versus AN where peak-peak cross-correlations are seen to decrease 
faster than both valley-valley and peak-valley cross-correlations. The enhanced cross-correlation for val- 
leys reflects the fact that, unlike resonant tunneling on peaks, elastic cotunneling relies on contributions 
from ~ E/A levels. In moving from one valley to the next, only one of the E/A levels is different, hence 
the similarity among neighbors. 

We thank I. Aleiner and L. Glazman for many valuable discussions. Work at Stanford support in 
part by tiie Office of Naval Research YIP program, the Army Research Office, the NSF-NYI program, 
and the A.P. Sloan Foundation. Work at UCSB supported by tiie AFOSR and QUEST. 
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Fig. 1 : (a,b) Conductance, G, as a function of magnetic field, B , and (a) as a function of gate voltage, 
Vg, and (b) for four particular paths. Paths used in (b) are shown as labeled white traces in (a). (c,d) 
Autocorrelations C, , (46) of conductance fluctuations for (c) all interpolated paths as a function of aver- 
age gate voltage along path, (Vg'j, and (d) for four traces, (d) Inset: Q ,(/46) near AB = 0. 
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Fig. 2: (a) Ensemble-averaged characteristic field B^. (solid) and average conductance (dashed) 

— — -2 

across peak-valley-peak for -14 independent data sets (dot 1), showing modulation of B^.. (b) B^ 
vs. {G)g for the same data. Diagonal line indicates B^^{G) = constant; horizontal line is saturation 
B^. = {k(f>f,/Ay^ = 0.003 mT"^ for Ej,< E. (c) Average B^ for three open dot configurations and peak 
and valley conductances (diamonds), from data in (a). Unaveraged B^. values for peaks (crosses) and 
valleys (circles) show spread in data. Open-dot 5^, is modified (by ~ 5-15 %) to reflect changes in dot area 
upon opening leads. Top inset: Schematic energy diagram of a blockaded dot. Bottom inset: micrograph 
of dot 1. 
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Fig. 3: Ensemble-averaged normalized conductance, g, (solid) for 81 statistically independent (a) peaks 
and (b) valleys as a function magnetic field for dot 3. Peak conductance gpg^k has a dip around B = 
with a width ~ associated with the breaking of time-reversal symmetry. No dip around = is seen 
in g valley- I^sets show the regions around B = 0. Unaveraged conductance (dashed) shows large fluc- 
tuations around average values. 
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Fig. 4: Average cross-correlation of conductance fluctuations, C(A5, AA^), for (a) peaks and (b) valleys, 
(c) Maximum cross-correlation at AB = 0, C(0, AA^) , for peak-peak (triangles), valley-valley (circles) and 
peak-valley (squares), showing peak-peak correlations decrease more quickly than the others, (d) Device 
parameters for the three dots: mean free path (1), Area (A) based on ~ 150 nm depletion around gates, 
mean level spacing ( A = 2;rh^ /m*A), charging energy {E(. = e^lC), and Thouless energy 
(£r=hVf/AV2). 
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